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We study effects of Kac-Baker long-range dispersive interaction (LRI) between particles on kink 
properties in the discrete sine-Gordon model. We show that the kink width increases indefinitely 
as the range of LRI grows only in the case of strong interparticle coupling. On the contrary, the 
kink becomes intrinsically localized if the coupling is under some critical value. Correspondingly, 
the Peierls-Nabarro barrier vanishes as the range of LRI increases for supercritical values of the 
coupling but remains finite for subcritical values. We demonstrate that LRI essentially transforms 
the internal dynamics of the kinks, specifically creating their internal localized and quasilocalized 
modes. We also show that moving kinks radiate plane waves due to break of the Lorentz invariance 
by LRI. 
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I. INTRODUCTION 



The effects of long-range dispersive interactions 
(LRFs) on the dynamics and thermodynamics of soliton- 
bearing systems have attracted a great deal of interest in 
the past decade Such attention is due to the fact 

that in realistic physical systems the interparticle forces 
are always long-ranged to some extent, and if the range 
of LRFs exceeds some critical value, they change soli- 
ton features qualitatively. In particular, the competition 
between short-range and long-range interactions in an- 
harmonic chains [0-^ and nonlinear Schrodinger (NFS) 
models brings into existence several types of the 

soliton states. In the nonlocal discrete NFS model two 
types of stable soliton states can coexist at the same 
excitation number p^ . In other words, there occurs a 
bistability phenomenon with a possibility of controlled 

Besides, the power law 
;braic soliton tails [E2pP] and 
can give rise to an energy gap between the spectra of 
plane waves and the soliton states . 

In the present paper we consider the effects of FRFs in 
discrete Klein-Gordon (KG) models. These models were 
successfully used in investigations of a number of physical 
phenomena such as dislocations in solids, charge-density 
waves, adsorbed layers of atoms, domain walls in fer- 
romagnets and ferroelectrics, crowdions in metals, and 
hydrogen-bonded molecules (see the review paper ]2^ ] 
for references). As it is known |16 2^ the interparticle 
interactions in many of these systems are substantially 
long-ranged. 



In the assumption of the harmonic interaction between 
particles the dimensionless Hamiltonian of the discrete 
KG model can be written in the form 



2 ^ ^ Jm.ni^ 



(1) 



switching between states 
FRI manifests itself in alge 



where Un is the displacement of the n-th particle from its 
equilibrium position and Jm.n is the coupling constant 
between particles n and m. 

As far as we know, until now there was only one inves- 
tigation jl^] of the KG model with the power law FRI. It 
was shown that the asymptotics of the kink shape as well 
as the interaction energy of the kinks are power law and, 
because of this, the dependence of the Peierls-Nabarro 
barrier versus the atom concentration is similar to the 
"devil's staircase" . 

But KG models with the exponential law (usually 
called Kac-Baker) FRI 



J(e" - l)e' 



■a rn— n 



(2) 



were believed to have been investigated in an exhaustive 
fashio n |l^ - [2^ ]. As early as 1981, Sarker and Krumhansl 
found ]17| an analytical kink solution for the 0* model. 
The width and the energy of the kinks were found to in- 
crease indefinitely as a decreases. An important role of 
the Kac-Baker FRI in thermodynamics of the system was 
also shown. Within the decade Woafo et al. considered 
in a series of papers p8|-p0[ the discreteness effects in the 



1 



same model. They have shown that the Peierls-Nabarro 
barrier vanishes as a — > 0. 

More recently the sine-Gordon (SG) model 



1 



cos Ur, 



(3) 



with Kac-Baker LRI (g) has been studied [£1 22 1 and all 
results of Ref. have been extended to this model. An 
implicit form for the kinks has been obtained and the 
kink energy and width have been found to grow to infin- 
ity as a ^ in Ref. |2|]. The thermodynamics of the 
system has been thereafter studied in Ref. ]2^. 

Thus, the investigations performed in Refs. |l7|~p^ give 
the impression that the Kac-Baker LRI always results, in 
the limit a — > 0, into infinite increasing of the kink width 
(and, therefore, vanishing of the Peierls-Nabarro barrier). 
However, closer inspection shows that this conclusion is 
proper for the case J > 0.5 only. 

In the present paper we explore the effects of the Kac- 
Baker LRI in the discrete sine-Gordon model fur- 
ther so that we could cover the case J < 0.5. What is 
more, we investigate the internal kink dynamics and the 
radiation of moving kinks. 

The paper is organized as follows. In Sec. |l|we derive 
the equation of motion of the system in the continuum 
limit using the technique of pseudo-differential operators. 
Then, in Sec. [II we solve this equation and obtain an 
implicit analytical form of the kink solution for arbitrary 
values of a and J. Turning back to the discrete case we 
calculate the form of the kinks numerically and compare 
it with the analytical solution. We show that in the case 
of J(e" -|- 1) < 1 the kinks are intrinsically localized. The 
calculation of the Peierls-Nabarro barrier as a function 
of a and J finishes the section. It turns out that the 
Peierls-Nabarro barrier vanishes in the limit a for 



J > 0.5 but remains finite for J < 0.5. In Sec. IV 
we develop a variational approach to the internal kink 
dynamics and demonstrate that LRI strongly enhances 
creation of kink's internal modes. Then we validate this 
result by direct numerical calculations. We show that 
similar to the non-sinusoidal Peyrard-Remoissenet poten- 
tial [^4|j25[ , the Kac-Baker LRI (||) with small a creates 
several kink's internal modes. By this means our results 
support the recent conclusion of Kivshar et al. that 
"the internal mode is a fundamental concept for many 
nonintegrable soliton models". Moreover, we show that 
for large values of J, for which kink's internal modes do 
not exist, the Kac-Baker LRI gives rise to pronounced 
quasilocalized modes inside of the phonon spectrum. In 
Sec. we show that due to break of the Lorentz invari- 
ance by LRI, there are no stationary moving kinks in the 
system; arbitrary moving kink will radiate plane waves 
with the wavelength proportional to its velocity. In Sec. 
[V^ we summarize and discuss the obtained results. 



II. EQUATIONS OF MOTION 

The Hamiltonian (|l|)-(||) generates the equation of mo- 
tion 



"diF 



ra^n 



sm Ur, 







(4) 



To obtain its solution analytically we pass to the con- 
tinuum limit treating n as a continuous variable n — > 
X = an, where a is the distance between particles. Thus, 
using 



Unit') u{x, t) 



(5) 



and keeping formally all terms in the Taylor expansion 
of ^ -vve can cast Eq. (^) in the operator form 



J(e° + l)smh^{ad^/2) 
sinh^(a/2) - sinh^(a9^/2) 



u + smu — , (6) 



where dx and dt are the derivatives with respect to x and 
respectively, and the identity 



— a\m\+arndx — 



_ cosh(a9x) — e " 



cosh(a) — cosh{adx) 



(7) 



has been used. 

In the approximation sinh(a92;/2) w adx/2 the equa- 
tion of motion (ra) takes on the form 



9? 



1 - (72^2 

where z = x/{a^) with 



sin u — 



v/J(e" + 1) 
2sinh(Q:/2) 



and 



J(e" + 1) 



(8) 



(9) 



Here the parameter ^ presents a measure for the soliton 
width — the continuum approximation should be good 
for large enough values of ^. 

Acting on Eq. (H) by the operator (1 — cr^d^) one can 
write the equation of motion in the differential form 



Utt 



smu 



o-'^Uzztt + cr^ (sin u)^ 



(10) 



which coincides with the equation derived in Ref. pl| . 
The authors of that paper found the form of the moving 
kink solution neglecting the term Uzztt- Thus they have 
found in fact an exact solution for immobile kinks but ap- 
proximate for moving ones. In Sec. we will show that 
the term Uzztt is responsible for the non-perturbative ra- 
diation of the kink. Indeed, due to break of the Lorentz 
invariance by LRI there are no stationary moving kinks 
in the system — moving kinks always radiate and even- 
tually stop. It is why in the next section we consider 
the immobile kinks only and write their exact shape in 
to some extent more simple and general form than that 
given in Ref. pi. 
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III. KINK'S STATIC PROPERTIES 

A. Analytical kink solution 

In this section wc obtain the immobile {dfu = 0) kink 
solution of Eq. (||). Denoting this solution as (j){y;B), 
where 



and 



z 2 . fa 
2/ = - = - smh — I a; 
a a V 2 



^2 = 1 + — = 1 + J(e" + 1) 



one can see that Eq. (|^) takes on the form 
(^^1 , / „N sin (b 



Q2 



1 



(11) 



(12) 



(13) 



Using a new variable v = sin(0/2) one can rewrite it in 
the polynomial form 



{l-V^){B^-2V^)Vyy 



(14) 



and multiplying by (B^ — 2ti^)/(l — w^)^ one can cast it in 
the form of the equation of motion of certain Hamiltonian 
system with the Hamiltonian 

Imposing the kink's boundary conditions {vy — > when 
w — > 0) we arrive at the constraint h = 0. Thus we obtain 
the equation 



«2(l-«2)(52 



(16) 



y (B2 - 2^2)2 

which after integration gives the kink solution of the form 



±(y-yo) = ylog(^^-^ 



log 



1 + ^l 



E 



m— 1 

where a new variable 

M = 



2m - 1 ' 



B^-v^ 



(17) 



(18) 



is introduced. Turning back to the function (/)(y; B) we 
obtain an exact form of the kink (positive sign in Eq. 
(|l7|)) or antikink (negative sign) centered at yo: 



(j){y; B) = 2 arcsin 



'1-^2^2 
1-//2 



(19) 



where the dependence of /i on y is determined by Eq. 
(0)- 

It can be checked that in the nearest-neighbors inter- 
action (NNI) limit (B — > oo) the solution reduces to the 
ordinary SG kink or antikink form 



u{z) = 4 arctan(exp(±(z — ^o))) 



(20) 



where we have used the variables u{z) and z defined in 
the previous section. 

Looking at Eq. ( p^ ) one can see that near the cen- 
ter of the kink, where /i is small, the first order term in 
the Taylor series vanishes at i? = \/2. In this case the 
derivative d^/dy = ±l/2/x^ goes to infinity in the center 
(/i 0) of the kink. It means that the slope of the kink 
becomes vertical for B = \/2 (see Fig. |l|). If B < \/2 the 
slope of the kink assumes negative values and the solu- 
tion ( pj| ) becomes 5-shaped (multivalued) and thus loses 
its physical meaning. 




-1 1 

2 -1/2 

z=(B -1) y 

FIG. 1. The kink shape predicted analytically for differ- 
ent values of the range of LRI: it is usual in the NNI limit 
(B = oo, full line), has the vertical slope in the critical case 
(J(e" -|- 1) = 1 or _B = \/2, dotted line), and is multivalued 
(5-shaped) in the supercritical cases of -B = 1.2 (dashed line) 
and _B = 1.1 (long-dashed line). 



Thus, returning to the initial physical parameters J 
and a we can state that there exists a critical value of J 
in the system: if J > 0.5 the value of B always exceeds 
\/2 and the form of the kink does not change drastically 
with a (see Fig. p[-a). It was this case which was studied 
in details in Refs. [ plp^ . But if J < 0.5 there is some 
critical value of a. for which B = \/2 and the transition 
from usual kinks to iS-kinks occurs when a decreases (see 
Fig. I^b). And now an interesting question should be 
raised: what is a physical single-valued analogue of the 
5-kink in the discrete case? 
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B. Numerical results 



The best remedy to answer the above question is to 
solve Eq. (|^) numerically. Since for a static solution Eq. 

turns itself into a system of N nonlinear algebraic 
equations (where A'^ is a number of particles), it is con- 
venient to use the Newton-Raphson iterations. To avoid 
perturbations due to boundary effects (we use a chain 
with fixed ends) the value of TV was chosen large enough 
(typically N = 500, but it has been extended to 1000 for 
broad kinks at small a and big J). The choice of the ini- 
tial kink form for zeroth iteration is not very important 
since for the given problem the Newton-Raphson itera- 
tions are very stable (but to be specific we used Eq. (|20| ) 
for this purpose). To obtain a stationary kink shape with 
the equilibrium positions of the particles in a chain 
we usually performed 7-20 iterations. 



eq 




-6 -4 -2 



them is excellent for all values of J. The only difference is 
that one should cut out the unphysical part of the 5-kink 
and replace it with a vertical slope to obtain the form of 
the kink in the discrete case (see Fig. ||). In our opin- 
ion this result can be understood by reference to a two- 
component kink structure. Indeed, as it was shown in 
Refs. [^U^ for the anharmonic chain with the Kac-Bakcr 
LRI between particles, the existence of two length scales 
results into two-component soliton structure, where the 
short-range component is dominant in the center of the 
soliton, while the long-range component (which can be 
properly described in the continuum approximation by 
Eq. (p^)) is dominant in the tails. 



eq 
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FIG. 3. The kink shape at small J and a: instead of the 
predicted analytically multivalued <S-shaped kink we obtain 
numerically the intrinsically localized kink. 
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FIG. 2. The comparison of the kink shape predicted an- 
alytically with that found numerically. Two cases must be 
distinguished: a) J > 0.5 for which the kink width increases 
indefinitely with decreasing of the a; b) J < 0.5 for which the 
kink width remains finite with decreasing of a. 

In Figs. we compare the form of the kinks found 
numerically for different J and a with the solution (p^- 
(H). One can see that at small a the agreement between 



^PN 




a 

FIG. 4. The dependence of the energy Epn of the 
Peierls-Nabarro barrier on the range a of LRI for different 
values of J. The barrier vanishes in the limit q — > for 
J > 0.5 but remains finite for J < 0.5. 

Thus, now we can conclude that in the discrete SG 
model the LRI affects the kinks in two opposite ways in 
relation to the value of J. When J > 0.5 the increas- 
ing of the range of LRI causes the increasing of the kink 
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width in agreement with the conclusion of Ref . . But 
for J < 0.5 the kinks become intrinsically localized as 
a — > 0. In the latter case the form of the kink is per- 
fectly described by the 5-kink in the tails with a vertical 
slope in the center. 

Numerical calculations show that in both cases the 
kink energy monotonically grows to infinity when a de- 
creases to zero. But the behavior of the Peierls-Nabarro 
barrier (defined as an energy difference of the kink cen- 
tered on a particle and the kink centered between parti- 
cles) completely correlates with the behavior of the kink 
form (see Fig. ^. When (for J > 0.5) the kink width 
grows in the limit a — > 0, the Peierls-Nabarro barrier 
vanishes. But when (for J < 0.5) the kink becomes 
intrinsically localized in this limit, the Peierls-Nabarro 
barrier remains finite. 



IV. KINK'S INTERNAL MODES 

In the previous section we were concentrating upon the 
static properties of the kinks. But a considerable interest 
is also attracted to the phonon spectrum affected by the 
presence of the kink. It is common knowledge that the 
influence of the kink is to some extent similar to that of 
the impurity affecting the phonon spectrum in a solid. 
Namely, not only quantitative changes in the spectrum, 
but qualitative ones consisting in emerging the localized 
modes with frequencies lying beyond the phonon band 
are to be expected and present special interest p5|-^ . 

There should be distinguished localized internal modes 
from the localized translational mode. The low-frequency 
translational mode (which is in the continuum limit the 
Goldstone mode associated with the translational invari- 
ance) is universally present in an arbitrary KG model. Its 
frequency is closely associated with the Peierls-Nabarro 
barrier considered above, so we shall not focus much at- 
tention on this mode in what follows; instead, we shall 
consider in detail the internal modes. These latter play 
an important role in the kink dynamics because they can 
temporarily store energy taken away from the kink's ki- 
netic energy, which can later be restored again in the 
kinetic energy. This gives rise to resonant structures in 
kink interactions ||3l|-p4|. 

The most extensively studied is the Rice internal mode 
prf which can be visualized as an oscillation of the kink's 
core width. Although this mode does not exist in the 
usual continuum SG model (which is integrable), even 
small perturbation of the model brings it into the exis- 
tence p5| , p6t . In particular, the Rice mode exists in the 
discrete SG model with NNI |25|. But the discreteness 
just changes the dispersion of the system and the dis- 
persive LRI under consideration affects it even greater. 
Thus, one might expect that the LRI will enhance the 
creation of the Rice internal mode in the system being 
considered. In the next subsection we recourse to a vari- 
ational approach and show that this is so indeed. Then 



we investigate kink's internal modes in more detail nu- 
merically. It turns out that at small a there can exist 
either several kink's internal modes below the phonon 
spectrum or pronounced quasilocalized modes inside the 
phonon spectrum. 



A. Variational approach 

When, as with the yj^-model, the Rice internal mode 
is pronounced it can be properly described by a varia- 
tional collective coordinates approach |^^. Proceeding 
from Eq. (^ and the identity 



^ ^ ^ ^ <-^m,n(Wm ^n) 



^ ^ ^ ^ <An,n(^m ^n) 



(21) 



one can pass to the continuum limit (||) and write the 
Hamiltonian (111) in the form 



H = 



1 r , {'^fduV ,^ 

„ . , / — TTTT / dy< - — + 1 — cos Ml 
2sinh(a/2) ''\2\dtJ ^ ' 



u{y,t) 



1 



92 



(22) 



where y = 2sinh(a/2) n. 

It should be emphasized that the Rice's collective co- 
ordinates approach p7[ cannot be used in our case. In- 
deed, choosing the trial function of the form (j){y/L{t); B) 
where (j){y;B) is the stationary kink solution (p^)-(p^) 
and L{t) is a time-dependent variational parameter (so- 
called "effective kink width" ) , we are unable to integrate 
analytically the long-range part of the Hamiltonian (|2^). 

Thus, we are forced to introduce another trial function. 
We call your attention to the fact (see Fig. |]) that the 
change of B in the kink solution <f>{y; B) changes a slope 
of the kink and its width as well. Therefore, to describe 
small-amplitude kink oscillations around its stationary 
form one can use equally well instead of Rice's trial func- 
tion a trial function of the form 



Un{t) = 0(2 sinh(a/2) n ; bit)) 



(23) 



where (j> is determined by Eqs. (|l7|)-([l9[) and b{t) is the 
time-dependent variational parameter. Then, using that 
(j){y; b{t)) is the solution of the equation 



1 - 92 

Uy 



sin 0(7/; 6) 
62 - 1 



(24) 



the integrals appearing in Eq. ( |22| ) can be taken analyt- 
ically: 

U{b)^ dy{l-coscb)=Ab-2log(^^y (25) 
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Vib) = - 

and 

Mib) = 



-J dy <!> sin 2b +{b'- I) logi^-^ 



(26) 



dy 



"1/6 

X / d/i- 





2bfi+{b^-l)\og[^ 



(1 - (62 - 2 + 62/i2) 



(27) 



where the last integral was taken analytically as well, 
yet its cumbersome expression prevents us from writing 
it down here in an explicit form. We just mention that 
M{b) grows to infinity as 6 ^ \/2 and has the asymp- 
totics 

(28) 

for big values of b. 

Thus the effective Hamiltonian of the system takes on 
the form 



H,ff ^ —^-—\l-M{b)(^) +W{b) 
2sinh(a/2) \2 \dt ) ^ ' 

where the potential energy of the kink 



(29) 



W{b) = U{b) 



62-1 



V{b) ~ Ek{B) 



;{b- Bf +0{b- Bf , 



(B2_l)2 

as it would be expected, is minimal in the point b = B 
■y/l + J(e" + 1), where it equals to the energy 

'B + V 



(30) 



EKiB) = 6B+ (B2 _3)iog 



B-1 



(31) 



of the stationary kink. The kink energy ( pi| ) was calcu- 
lated first (although in a more bulky formjin Ref. pl[|. 

For small deviations of 6 from B (in the harmonic ap- 
proximation) the effective Hamiltonian (|2^) generates the 
equation of motion 



— + n^]{bit)-B) = o 



(32) 



being the equation of motion for the harmonic oscillator 
with the frequency 



1/2 



M{B){B^ ~ 1)2_ 
whose dependence on B is depicted on Fig. ||. 



(33) 
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variational approach 



6 o ❖ numerical (for J=0.6) 



1 



3 
B 



FIG. 5. The dependence of the frequency Q. of the 
Rice's kink internal m ode on the parameter of the nonlocality 
B = yTT"j(e^TTy: predicted analytically by Eq. || (full 
line) and found numerically for J = 0.6 (diamonds). 

Thus, a slightly excited kink will oscillate around its 
stationary shape with the frequency fi, depending upon 
the parameter B which specifies the nonlocality of the 
system. When B decreases (nonlocality grows) the fre- 
quency ri also decreases (57 — > when B V2). Al- 
though Eq. ( p3| ) is not fully precise even in the contin- 
uum limit (to improve it an interaction with phonons 
should be taken into account |2^) it is in a good agree- 
ment with the numerical calculations (see Fig. ^ which 
are described below. 



B. Numerical results 



The exposed variational approach (even in a form com- 
plicated by introducing of several time-dependent param- 
eters) permits one to investigate only a limited number 
of oscillatory modes. To overview rather the entity of the 
whole phonon spectrum one should deal with the initial 
set of the equations of motion (^). Specifically, when all 
the equilibrium positions u^j'^ of the particles in a chain 
with a kin k become known by the method described in 
Sec. IIIB , we can study the spectrum of small-amplitude 
oscillations around this state by looking for a solution of 
Eq. (0) in the form 



(34) 



We assume that the deviations w„ of the particles from 
the kink shape are sufficiently small (d„ <C a) and ignore 
all nonlinear terms in the equations of motion for Vn ■ In- 
troducing this ansatz into Eq. (^) , we obtain a system of 
linear equations that can be written in a matrix form 



Dv = n^v 



(35) 
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where v = {vn} and a symmetric matrix D is the dynam- 
ical matrix of the lattice in the presence of a kink with 
the components: 

Dn,n= 2J + cos<'i , 

By virtue of the fact that the matrix D is symmetric all 
its eigenvalues are real. They give us the frequencies of 
the small-amplitude oscillations around the kink while 
the corresponding wavevectors describe the spatial pro- 
file of each mode. The eigenvalues and eigenvectors of the 
matrix D were calculated numerically using the House- 
holder matrix reduction to tridiagonal form and the QL 
diagonalization algorithm The obtained results are 
presented in Figs. ^ - In the following discussion 
of these results, we assume that the reader is familiar 
with Ref. 1^ where the basic features of the kink's lin- 
ear spectrum are comprehensively expounded. 




In Fig. ^ we compare the kink's linear spectrum ri( J) 
for a = 0.5 with the spectrum obtained in Ref. for 
the NNI limit (a = oo). As it has there been shown, in 
the discrete SG model with the interaction between near- 
est neighbors aside from the low-frequency translational 
localized mode there exist (slightly expressed in the in- 
terval 0.27 ^ J ^ 1.59) the Rice internal mode. One 
can see from Fig. ^ that for a = 0.5 the situation re- 
mains qualitatively the same, except that the Rice mode 
becomes evidently pronounced. Thus, in agreement with 
the results of the variational approach, the LRI strongly 
enhances creation of the Rice mode (see also Fig. ||). 




FIG. 7. Spectrum of small-amplitude excitations around a 
kink as a function of the nonlocality parameter a for J — 0.6. 




J 

FIG. 6. Spectrum of small-amplitude excitations around 
a kink as a function of the coupling parameter J for a = oo 
and a — 0.5. 



Furthermore, when a yet decreases, there occur new 
qualitative phenomena. Foremost, additional localized 
internal modes are split out from the bottom of the 
phonon spectrum. This process is easily observable in 
Fig. ^ where the kink's linear spectrum Q as a, func- 
tion of a is plotted for J = 0.6. One can see that the 
number of localized modes grows indefinitely as a van- 
ishes. In particular, there are four localized modes at 
a = 0.2 (see Fig. ||-a) and seven at a = 0.1 (see Fig. ||- 
b). All the localized internal modes are best pronounced 
around J ~ 0.6. Closer analytical examination of Eqs. 
(^)-(|36|) shows that all eigenstates are either symmetric 
or antisymmetric. One can state on the strength of the 
numerical calculations that the symmetric and antisym- 
metric states are always alternating, starting with the 
symmetric translational mode (Mode 1 in Fig. ^ and 
antisymmetric Rice mode (Mode 2 in Fig. |9|). 
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FIG. 8. Spectrum of small-amplitude excitations around a 
kink as a function of the coupling parameter J for: a) a = 0.2 
and b) Q = 0.1. 



Another interesting feature that appears in the kink's 
linear spectrum at small a is that there exist an appar- 
ent transformation of the phonon spectrum (in the plane 

— J) along an imaginary line (let us call it "ridge") 
which extends the localized internal modes inside the 
phonon spectrum in the direction of larger J (see Fig. 
g). This ridge is evidently parallel to the upper cut- 
off frequency of the phonon band, with the spectrum 
curves looking like essentially different left or right from 
the ridge. In the domain right to the ridge each sym- 
metric eigenstate is equidistant from both nearest anti- 
symmetric eigenstates. In the domain left to the ridge 
each symmetric eigenstate is paired with an appropri- 
ate antisymmetric eigenstate: their frequencies become 
almost coincident. The shapes of the symmetric and an- 
tisymmetric states in a pair are also in a close correla- 
tion: ~ 9{n — no)v^ , where no is the position of 
the kink center and d{x) is the Heaviside function. The 
ridge is shifted towards the upper cut-off frequency as 
a grows and completely disappears (or rather coincides 
with the upper edge of the phonon band) in the NNI limit 
(a = oo). 
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FIG. 9. Example of the shape of the lowest localized 
modes: the translational mode (Mode 1), the Rice mode 
(Mode 2), and so on. 
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FIG. 10. The shape of the delocalized and quasilocalized 
states inside the phonon band for J = 1.25 and a = 0.1. 

The most interesting phenomenon however concerns 
the shape of the eigenstates in the vicinity of the ridge 
(see Fig. Ml One can see that both above (e.g., mode 
40 in Fig. |10|) and below (e.g., mode 8 in Fig. ^ the 
ridge the eigenstates are delocalized as it is expected to 
be in the continuum spectrum. But the shape of the 
eigenstates undergoes essential changes when approach- 
ing the ridge frequency. Namely, such eigenstates prac- 
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tically have no distortion in the oscillatory subspace and 
thus are quasilocalized (see modes 11-13 in Fig. |l^). It 
should be indicated a similarity between these quasilocal- 
ized states and the "exotic states" investigated in Refs. 
p6|-^, although wc could not yet provide an in-depth 
analysis of this analogy. 
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FIG. 11. The density of the eigenstates in the phonon band 
as a funcion of the eigenfrequency for J = 2.0 and a = 0.1. 

It is notable that the quasilocalization of the eigen- 
states in the vicinity of the ridge frequency is accompa- 
nied by pronounced densening of the eigenstates which 
can be viewed in Fig. ^ where we plot (at fixed a and 
J) the density of the eigenstates as function of frequency. 
Because of this peculiarity inside the phonon spectrum, 
the ridge of quasilocalized states should be observable ex- 
perimentally. We expect that similar to the localized in- 
ternal modes, the quasilocalized states could play a part 
in the kink dynamics. 



V. RADIATION OF MOVING KINKS 

All the preceding was devoted to the properties of im- 
mobile kinks. In this section we go out of this restric- 
tion and consider effects of the LRI on the moving kinks. 
To avoid discreteness effects we assume that J > 0.5 
and ^ ^ 1, so that the kinks are broad and the Peierls- 
Nabarro barrier can be neglected. In other words we as- 
sume that the kinks are properly described by Eq. (||) . It 
is well known that in the limit of short-ranged dispersion 
(cr = in Eq. (|^)) the kinks can move without changing 
their shape and velocity. In fact, it is a consequence of 
the Lorentz invariance of Eq. (||) for a — Q. But this 
invariance is broken by the LRI (as well as it would be 
broken by the discreteness). Therefore, one might expect 
that similar to the discreteness the LRI causes a ra- 
diation of moving kinks. Indeed, this effect was recently 
shown to exist for fluxons in a superlattice of Josephson 
junctions pO| which is described by another variant of 
nonlocal SG equation. In what follows we extend the 



results of Ref. [Q to our model. 

To investigate propagating solutions we introduce the 
transformation to the moving frame of reference in which 
the center of the kink is at rest: 



(t — vz) 7(f) 



(37) 



where j{v) = l/\n~~v^ is the Lorentz factor and v is 
the dimcnsionlcss kink velocity. Applying this transfor- 
mation to Eq. (^) we obtain 



-f\dl-2vdrdc+v^dl)u{C,T) 



1 + a'^^'^k'^iv) 
where the notation 

-.2 



k\v)^2vdrdQ-dl-v^dl 



was used. We introduce 

"(C,^) = 0(CM;S) 



/(C,r) 



(38) 



(39) 



(40) 



where 6{(/aj;B) is the static solution given by Eqs. 
(|l7|)-([l9|) and the function /(C,t) describes the change 
of the kink shape and the radiation. Inserting Eq. ( |40| ) 
into Eq. (^8|) in the linear approximation with respect to 
/(C,t) we obtain for the Fourier transform of f{C,T) 



Kk,T) 



1 

2^ 



-^■';/(C,T)dC 



(41) 



the inhomogeneous differential equation 

^\dl-2ivkdr)f + n^{k,v)f 



1 + 720-2^2 



-(1 



cos ( 



1/ 



1 + -f'^a'^k'^ , 



I 



where 



k^{v)^k^ + 2ivkdr-v^dl 



and the function 



n'^{k,v) = l-72t)2A;2 + 



72fc2 



(42) 



(43) 



(44) 



1 + "f^a^k^ 

is the dispersion relation in the moving frame of refer- 



ence. The right-hand-side of Eq. (42) is the source of 
the radiation. It vanishes when the kink velocity v = 0. 

The equation Q{k, v) = has two real roots k = ±kr. 
The existence of these roots means ^,^,|l^ that waves 
with wavenumbers ± kr will be resonantly excited by a 
kink, forming an oscillatory tail. We restrict ourselves 
to the investigation of small kink velocities for which the 
resonant wavenumber kr can be expressed as follows 



^ V 



(45) 
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In the resonant region A; ~ ±kr we can neglect the sec- 
ond derivative f because as it will be seen later the 
function / slowly changes with r. We use the expansion 



(fc, v) ~ 2"f V kr [kr T k) for k ~ ±kr 



(46) 



and replace the right-hand side of Eq. ( [42| ) by its value 
at fc = ±fcr- We represent also the function f{k,T) as a 
sum 



f^f+ + f- , 



(47) 



where the function /+ (/-) differs from zero only in the 
close vicinity of the wave vector kr {—kr). Taking into 
account Eqs. (|46|)-(|^ we obtain from Eq. (^) that the 
functions f± satisfy the equation 



idr f±+ v{k =F kr)f± 



where 



A = (dcmr) 



sech 



ttB 
27 



(48) 



(49) 



Returning to the real space t) we obtain from Eq. 

(H) 



{dr - vd^)f± ± i 



f± = ^S{0. (50) 



VI. CONCLUSIONS 

We have investigated the effects of the Kac-Baker long- 
range interaction on the kink's properties in the discrete 
sine-Gordon model. We have obtained an implicit form 
for the kink's shape and energy and have shown that the 
kink width increases indefinitely as a vanishes only in 
the case of strong interparticle couphng (J > Jcr = 0.5). 
On the contrary, the kink becomes intrinsically localized 
for J(e" -I- 1) < 1. Accordingly, we have shown that the 
Peierls-Nabarro barrier vanishes as a ^ for supercrit- 
ical values of the coupling J but remains finite for sub- 
critical values. We have developed a new variant of the 
collective coordinates variational approach for investiga- 
tion of the internal kink's dynamics and have shown that 
the Kac-Baker LRI essentially enhances creation of the 
kink's internal localized modes. We have demonstrated 
numerically that indefinite number of the localized in- 
ternal modes came into existence when a approaches 
zero. We have revealed an existence of the quasilocal- 
ized states inside the phonon spectrum for small a and 
large J. We have described their properties, in particu- 
lar, a pronounced densening of the phonon spectrum in 
the vicinity of the quasilocalized states. We expect that 
similar to the localized internal modes, the quasilocal- 
ized states could play a part in the kink dynamics. We 
have investigated moving kinks in the continuum limit 
and have shown that, due to break of the Lorentz invari- 
ance by the LRI, they always radiate plane waves with 
the wave length proportional to the kink velocity. 



With the initial condition f{(,T)/r=o = 0, the solution 
of Eq. (M) becomes 



f = A {eic + vt) - eiO) cos {kr c - x(C)) , (5i) 



where 



f\^- cos cj,{C/ar,B)] dC 
2Bv Jo 



Bv 



tanh(C) 



(52) 
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is the phase function. It can be seen that \kr(\ ^ xiC) 
when ICI > 1. Therefore, we can neglect the phase func- 
tion x{0 returning to the original variables (z, t) we 
obtain that the radiation of the kink which moves with 
a small velocity v is given by the function 

f{z, t)^ A [e{z) - e{z - vt)] cos[{B/v) z-Bt], (53) 

with the amplitude of the radiation A which exponen- 
tially decreases when the kink velocity tends to zero. 

Thus, now one can conclude that the Kac-Baker LRI 
produces non-perturbative radiation with the wave length 
A — 2t: /kr ~ V that emerges in the rear of the kink. 
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